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Abstract 

We consider the problem of predicting a real random variable from a func- 
tional explanatory variable. The problem is attacked by mean of nonpara- 
metric kernel approach which has been recently adapted to this functional 
context. We derive theoretical results by giving a deep asymptotic study of 
the behaviour of the estimate, including mean squared convergence (with rates 
and precise evaluation of the constant terms) as well as asymptotic distribu- 
tion. Practical use of these results are relying on the ability to estimate these 
constants. Some perspectives in this direction are discussed. In particular a 
functional version of wild bootstrapping ideas is proposed and used both on 
simulated and real functional datasets. 
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1 Introduction 



Functional data are more and more frequently involved in statistical problems. De- 
velopping statistical methods in this special framework has been popularized during 
the last few years, particularly with the monograph by Ramsay & Silverman (2005). 
More recently, new developments have been carried out in order to propose non- 
parametric statistical methods for dealing with such functional data (see Ferraty 
& Vieu, 2006, for large discussion and references). These methods are also called 
doubly infinite dimensional (see Ferraty & Vieu, 2003). Indeed these methods deal 
with infinite-dimensional (i.e. functional) data and with a statistical model which 
depends on an infinite-dimensional unknown object (i.e. a nonparametric model). 
This double infinite framework motivates the appellation of Nonparametric Func- 
tional Statistics for such kind of methods. Our paper is centered on the functional 
regression model : 

Y = r{X) + error, (1) 

where F is a real random variable. A:" is a functional random variable (that is, X 
takes values in some possibly infinite-dimensional space) and where the statistical 
model assumes only smoothness restriction on the functional operator r. At this 
point, it worth noting that the operator r is not constrained to be linear. This is a 
Functional Nonparametric Regression model (see Section|2for deeper presentation). 

The aim of this paper is to extend in several directions the current knowledges 
about functional nonparametric regression estimates presented in Section |21 In 
Section 13.11 we give asymptotic mean squared expansions, while in Section 13.21 the 
limiting distribution is derived. The main novelty/difficuly along the statement of 
these results relies on the exact calculation of the leading terms in the asymptotic 
expressions. Section HI points out how such results can be used when the functional 
variable belongs to standard families of continuous time process. The accuracy of 
our asymptotic results leads to interesting perspectives from a practical point of 
view: minimizing mean squared errors can govern automatic bandwidth selection 
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procedure while the hmiting distribution of the error is a useful tool for building 
confidence bands. To this end, we propose in Section El a functional version of the 
wild bootstrap procedure, and we use it, both on simulated and on real functional 
datasets, to get some automatic rule for choosing the bandwidth. The concluding 
section IHl contains some important open questions which emerge naturally from the 
theoretical results given in this paper, such as the theoretical study of the accuracy 
of the functional wild bootstrap procedure used in our applications. 

2 Kernel nonparametric functional regression 

2.1 The model 

The model is defined in the following way. Assume that [Yi, Afj) is a sample of n 
i.i.d. pairs of random variables. The random variables Yi are real and the A^'s are 
random elements with values in a functional space S. In all the sequel we will take 
for S a separable Banach space endowed with a norm ||-||. This setting is quite 
general since it contains the space of continuous functions, spaces as well as more 
comphcated spaces like Sobolev or Besov spaces. Separability avoids measurability 
problems for the random variables A'j's. The model is classically written : 

Yi = r{Xi)+ei i = l,...,n, 

where r is the regression function mapping S onto M and the ^j's are such that for 
all i, E{ei\Xi) = and E{e^\Xi) = al{X) < oo. 

2.2 The estimate 

Estimating r is a crucial issue in particular for predicting the value of the response 
given a new explanatory functional variable Xn+i- However, it is also a very del- 
icate task because r is a nonlinear operator (from S into M) for which functional 
linear statistical methods were not planned. To provide a consistent procedure to 
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estimate the nonlinear regression operator r, we propose to adapt the classical finite 
dimensional Nadaraya-Watson estimate to our functional model. We set 

^ ELiY>^K {h-^ - x\\) 

El=iK{h'^\\x,-x\\) ■ 

Several asymptotic properties of this estimate were obtained recently. It turns out 
that the existing literature adresses either the statement of upper bounds of the 
rates of convergence without specification of the exact constants (see Chapter 6 in 
Ferraty & Vieu, 2006), or abstract expressions of these constants which are unusable 
in practice (as for instance in the recent work by Masry, 2005, which has been 
published during the reviewing process of this paper). Our aim in this paper is to 
give bias, variance, means square errors and asymptotic distribution of the functional 
kernel regression estimate with exact computation of all the constants (see Section 
Ej). We will focus on practical purposes in Sectional 

Several assumptions will be made later on the kernel K and on the bandwidth 
h. Remind that in a finite-dimensional setting pointwise mean squared error (at x) 
of the estimate depends on the evaluation of the density (at x) w.r.t. Lebesgue's 
measure and on the derivatives of this density. We refer to Schuster (1972) for an 
historical result about this topic. On infinite-dimensional spaces, there is no measure 
universally accepted (as the Lebesgue one in the finite-dimensional case) and there 
is need for developping a "free-density" approach. As discussed along Section HI the 
problem of introducing a density for X is shifted to considerations on the measure 
of small balls with respect to the probability of X. 

2.3 Assumptions and notations 

Only pointwise convergence will be considered in the forthcoming theoretical results. 
In all the following, x is a fixed element of the functional space S. Let (p be the real 
valued function defined as 

ip{s) = E[{r{X)-rix))\\\X-x\\=s], 
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and F be the c.d.f. of the random variable \\X — x\\'- 

F{t)^P{\\X-x\\<t). 

Note that the crucial functions ip and F depends implicitely on x- Consequently we 
should rather note them by ip^ and F^ but, as x is fixed, we drop this index once 
and for all. Similarly, we will use in the remaining the notation cr^ instead of cr^(A'). 
Let us consider now the following assumptions. 

HO : r and are continuous in a neighborhood of x, and F(0) = 0. 

HI : v7'(0) exists. 

H2 : The bandwidth h satisfies lim„^oo/i = and lim„^oo '"-^(^) = oo, while 
the kernel K is supported on [0, 1], has a continuous derivative on [0, 1), K'{s) < 
and X(l) > 0. 

Assumptions HO and H2 are clearly unrestrictive, since they are the same as 
those classically used in the finite-dimensional setting. Much more should be said 
on assumption HI. Note first that, obviously, p (0) = 0. It is worth noting that, 
whereas we could expect assumptions on the local regularity of r (as in the finite- 
dimensional case), hypothesis HI skips over that point and avoids to go into formal 
considerations on differential calculus on Banach spaces. To fix the ideas, if we 
assumed differentiability of r, we would get by Taylor's expansion that 

t{X)-t (x) = (r' {x),X-x)+o{\\X- xll) , 

where r' (x) G being the conjugate space of S and (•, •) being the duality 

bracket between £ and 8* . In this context, a non trivial link would appear between 
</?'(0) and r'(x) through the following relation: 

lim E 

s-»0 

Indeed, even if the link between the existency of r'(x) and of (^'(0) is strong, one 



r' (x) , 
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can build counter-examples for their non equivalence (these counter-examples are 
available on request but they are out of the main scope of this paper). In the 
perspective of estimating the constants given in Theorem ^ it will be easier to 
estimate v^'(O) (for instance by using r) than the operator r'(x). Therefore, we prefer 
to express computations by mean of v''(0) instead of r'{x)- This has the additional 
advantage to produce more readable writings. Consequently, the differentiability of 
r is not needed. 

Let us now introduce the function Th defined for all s G [0, 1] as: 

Tk (s) = = P{\\X~ x\\ < hs\ \\X - x\\ < h) , 

for which the following assumption is made: 

H3 : For all s G [0, 1] , Th (s) ^ tq (s) as /i ^ 0. 
Note that the function Th is increasing for all h. The measurable (as the pointwise 
limit of the sequence of measurable functions Th) mapping tq is non decreasing. Let 
us finally mention that this function tq will play a key role in our methodology, in 
particular when we will have to compute the exact constant terms involved in our 
asymptotic expansions. For the sake of clarity, the following proposition (whose a 
short proof will be given in the Appendix) will explicit the function tq for various 
cases. By l]o,i](-) we denote the indicator function on the set ]0, 1] and (■) stands 
for the Dirac mass at 1. 

Proposition 1 

i) If F (s) ~ Cs'^ for some 7 > then ro(s) = s'^ . 

a) If F (s) ~ Cs'^ |lns|'' with 7 > and k > then tq (s) = s'^. 

in) If F (s) ~ Cis''' exp (— C2/s^) for some p > and some 7 > then tq (s) = 

Si (s) . 

iv) If F (s) ~ C/ |lns| then Tq (s) = l]o,i](s). 

A deeper discussion linking the above behavior of F with small ball probabilities 
notions will be given in Section |3] 
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3 Asymptotic study 



In both following subsections we will state some asymptotic properties (respectively 
mean squared asymptotic evaluation and asymptotic normality) for the functional 
kernel regression estimate f. 

It is worth noting that all the results below can be seen as extensions to functional 
data of several ones already existing in the finite-dimensional case (the literature is 
quite extensive in this field and the reader will find in Sarda & Vieu (2000) deep 
results as well as a large scope of references). With other words, our technique 
for proving both Theorem Q and Theorem |21 is also adapted to the scalar or vector 
regression model since the abstract space £ can be of finite dimension (even, of 
course, if our main goal is to treat infinite-dimensional cases). Moreover, it turns 
out that the transposition to finite-dimensional situations of our key conditions (see 
discussion in Section |3] below) becomes (in some sense) less restrictive than what is 
usually assumed. With other words, the result of Theorem^ and Theorem |^ can be 
directly applied to finite-dimensional settings, and will extend the results existing 
in this field (see again Sarda & Vieu, 2000) to situation when the density of the 
corresponding scalar or multivariate variable does not exist or has all its successive 
derivatives vanishing at point x (see discussion in Section lOj) . 

All along this section we assume that assumptions H0-H3 hold. Let us first 
introduce the following notations: 
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3.1 Mean Squared Convergence 

The following result gives asymptotic evaluation of the mean squared errors of our 
estimate. The asymptotic mean squared errors have a standard convex shape, with 
large bias when the bandwidth h increases and large variance when h decays to zero. 
We refer to the Appendix for the proof of Theorem 

Theorem 1 When H0-H3 hold, we have the following asymptotic developments : 
Erix) -r{x) = ^' (0) ^h^O {{nF (h))-') +o{h), (2) 

and 

3.2 Asymptotic Normality 

Let us denote the leading bias term by: 

= ^' (0) ^^h. 

Before giving the asymptotic normality, one has to be sure that the leading bias 
term does not vanish. This is the reason why we introduce the following additional 
assumption: 

H4 : (^'(0) 7^ and Mq > 0. 
The first part of assumption H4 is very close to what is assumed in standard finite- 
dimensional literature. It forces the nonlinear operator r not to be too smooth (for 
instance, if r is Lipschitz of order /3 > 1, then v''(0) = 0). The second part of as- 
sumption H4 is specific to the infinite-dimensional setting, and the next Proposition 
121 will show that this condition is general enough to be satisfied in some standard 
situations. This proposition will be proved in the appendix. 

Proposition 2 

i) If tq{s) 7^ l]o,i](s) and tq is continuously differentiahle on (0, 1), then Mq > for 
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any kernel K satisfying H2. 

a) If tq{s) = Si{s), then Mq > for any kernel K satisfying H2. 

To emphasize the interest of these results, they should be combined with those 
of Proposition H Note that the result i) includes the well-known family of processes 
for which tq{s) = s'^ (see Proposition ^i), that is those whose distributions admit 
fractal dimensions (see Section l^?^ . The second case when tq{s) = Si{s) (for which 
a particular case is given in Proposition ^zii) corresponds to nonsmooth processes 
(see Section f4.H) . These two cases cover a large number of situations. However, 
if a more general function tq has to be used, one can make additional hypotheses 
on the kernel K. In particular, if K is such that Vs G [0,1], {sK{s))' > and 
ro(s) 7^ l]o,i](s), then Mq > 0, which covers the case of the uniform kernel. More 
complicated kernel functions K would lead to more technical assumptions linking 
K with tq. It is out of purpose to give these tedious details (available on request) 
but let us just note that the key restriction is the condition ro(s) 7^ l]o,i](s) (else we 
have Mo = 0). 

Moreover, since the rate of convergence depends on the function F (h) and for 
producing a reasonably usable asymptotic distribution it is worth having some esti- 
mate of this function. The most natural is its empirical counterpart: 



The pointwise asymptotic gaussian distribution for the functional nonparametric 
regression estimate is given in Theorem|21below which will be proved in the appendix. 
Note that the symbol ^ stands for "convergence in distribution". 

Theorem 2 When HO-H4 hold, we have 



A simpler version of this result is stated in Corollary Q below whose proof is obvious. 
The key-idea relies in introducing the following additional assumption: 



F{h) 



i^{^■■\\X,-x\\<h) 



n 
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H5 : limn^oo h>ynF{h) = 
which allows to cancel the bias term. 

Corollary 1 When H0-H5 hold, we have 

\lnF{h) {r{x) -r{x)) M (0, 1) . 

In practice, the constants involved in Corollary Q need to be estimated. In order 
to compute explicitely both constants Mi and M2, one may consider the simple 
uniform kernel and get easily the following result: 

Corollary 2 Under assumptions of Corrollary^ if K{.) = l[o_i](.) and if d1 is a 
consistent estimator of a^, then we have: 



^2 



There are many possibilities for constructing a consistent conditional variance esti- 
mate. One among all the possibilities consists in writing that 

al{x) = E{{Y-r{X)f\X = x), 

= E{Y^\X = x)-{E{Y\X = x))\ 

and, by estimating each conditional expectation with the functional kernel regression 
technique. 



4 Some Examples of small ball probabilities 

The distribution function F plays a prominent role in our methodology. This appears 
clearly in our conditions (through the function tq) and in the rates of convergence 
of our estimate (through the asymptotic behavior of the quantity nF{h)). More 
precisely, the behaviour of F around turns out to be of first importance. In 
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other words, the small ball probabilities of the underlying functional variable X 
will be determining. In order to illustrate our ideas and to connect with existing 
probabilistic knowledges in this field, let us now just discuss how F (and hence Tq) 
behave for different usual examples of processes X valued in an infinite-dimensional 
space. 

4.1 Nonsmooth processes 

Calculation of the quantity Pdl^Y — x|| < s) for "small" s (i.e. for s tending to 
zero) and for a fixed x is known as a "small ball problem" in probability theory. 
This problem is unfortunately solved for very few random variables (or processes) 
X, even when x = 0. In certain functional spaces, taking x 7^ yield considerable 
difficulties that may not be overcome. Authors usually focus on gaussian random 
elements. We refer to Li & Shao (2001) for a survey on the main results on small 
ball probability. If is a gaussian random element on the separable Banach space 
£ and if x belongs to the reproducing kernel Hilbert space associated with X, then 
the following well-known result holds: 



So, the small ball problem at any point x may be shifted to a small ball problem 
at 0. Moreover, can be precised in a few situations. For instance, Mayer-Wolf & 
Zeitouni (1993) investigate the case when A" is a one- dimensional diffusion process 
and X satisfies some conditions (see Mayer- Wolf & Zeitouni, 1993, pi 5). They also 
briefly mention the non gaussian case (see Mayer- Wolf & Zeitouni, 1993, Remark 3, 
pl9) but many other authors have considered different settings (see Ferraty et ai, 
2005, for a large discussion and references therein). As far as we know, the results 
which are available in the literature are basically all of the form: 



P{\\X-x\\<s)r.C^P{\\X\\<s) 



as s 



0. 



(4) 



Pi\\X-x\\<s)r^C, 



x 



s-" exp (-C/s^) 



(5) 
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where a,/3, c,^ and C are positive constants and ||-|| may be a sup, a W , a Besov 
norm ... The next remark is a direct consequence of Proposition^ It proves that 
non-smooth processes may satisfy the assumptions needed to get the asymptotic 
expansions of previous sections. 

Remark 1 In the case of "non-smooth" processes defined by (0) we have tq (s) = 
5i(s). In addition, condition nF{h) +00 (in H2) is checked as soon as h = 
A/ (logn)^^^ for A large enough. 

4.2 Fractal (or geometric) processes 

Another family of infinite dimensional processes is the class of fractal processes for 
which the small ball probabilities are of the form 

P{\\X-x\\<s)^c'^s\ (6) 

where and 7 are once again positive constants. Like above, it is elementary to 
get the following result from Proposition^ 

Remark 2 Under we have to(s) = s"^ while the condition nF{h) +00 (in 
H2) is satisfied as soon as h = An~^ for B small enough. 

4.3 Back to the finite dimensional setting 

Finally, it is important to note that a special case of fractal processs is given by the 
usual multivariate case (that is, by the case when £ = W). The following result is 
obvious for the uniform norm on MP and extends directly to any norm, since all of 
them are equivalent in finite dimension. 

Remark 3 If S = MP , then any random variable X on MP which has a finite and 
non zero density function at point x satisfies (0j with 7 = p. 

From Remarks |21 and |21 it is clear that all the results of Section |21 apply in a 
finite dimensional setting. Besides, the assumptions needed for Theorems Q and 
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121 are weaker than those described in Remark |21 since there is no need to assume 
the existence of a density for X. In this sense, our results extend the standard 
multivariate literature (see the discussion at the beginning of Section |2)). 

5 Perspectives on bandwidth choice 
5.1 Introduction 

The asymptotic results presented in the previous Section |21 are particularly appeal- 
ing because, in addition to the specification of the rates of convergence, the exact 
constants involved in the leading terms of each result are precised. This is particu- 
larly interesting in practice. Let us focus now on the mean squared errors expansion 
given in Section 13.11 In fact. Theorem [T] could give clues for possible automatic 
bandwidth choice balancing the trade-off between variance and squared-bias effects. 
However, the constants are unknown in practice which could seem to be a serious 
drawback for practical purposes. This general problem is well-known in classical 
nonparametric statistics, but in our functional context this question gets even more 
intricate because of the rather complicated expression of Mq, Mi and M2. An ap- 
pealing way to attack the problem is to use bootstrap ideas. In Section 15.21 we 
propose a track for building a functional version of the so-called wild bootstrap. 
We will show in Section through some simulated examples, how this functional 
wild bootstrap procedure works on finite sample sizes for choosing automatically an 
optimal bandwidth. A case study, based on spectrometric functional data coming 
from the food industry, will be shortly presented in Section 15.41 

At this stage it is worth noting that we have no asymptotic support for this 
functional bootstrapping procedure. This open question will be one of the main 
point discussed in the concluding Sectional 
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5.2 A Functional version of the wild bootstrap 

Basically, when using bootstrapping techniques one expects to approximate directly 
the distribution of the error of estimation without having to estimate the lead- 
ing terms involved in some asymptotic expansion of this error. In standard finite- 
dimensional problems (that is, when the variable X is valued in MP) , a so-called wild 
bootstrap has been constructed for approximating the distribution of the error of 
estimation in kernel nonparametric regression. We refer to Hardle (1989) and Hardle 
& Marron (1991) for a previous presentation of the wild bootstrap in nonparameric 
regression. A selected set of additional references would include Mammen (2000) 
for the description of the state of art on nonparametric bootstrapping, Mammen 
(1993) for a large study of wild bootstrap, and Hardle, Huet & Jolivet (1995) for 
specific advances on wild bootstrap in (finite-dimensional) nonparametric regression 
setting. 

The main interest of this kind of bootstrap relies on a resampling procedure of 
the residuals which makes it easily adaptable to our functional setting. Precisely, an 
adaptation to our functional setting could be the following functional wild bootstrap 
procedure: 

i) Given the estimate r constructed with a bandwidth h. compute the residu- 
als ei — Hi — r{Xi), and construct a sequence of bootstrapped residuals such 
that each e* is drawn from a distribution G* which is the sum of two Dirac 
distributions : 



Such a distribution ensures that the first three moments of the bootstrapped 
residuals are respectively 0, and (see Hardle & Marron, 1991, for details). 

ii) Given the bootstrapped residuals e*, and using a new kernel estimate f which 
is defined as f but by using another bandwidth g, construct a bootstrapped 
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sample {X*,Y*) by putting 



X* = Xi and Y; = r{Xi) + e*. 



iii) Given the bootstrapped sample {X*,Y*), compute the kernel estimate r*(x) 
which is defined as f (with the same bandwidth h) but using the bootstrapped 
sample {X*,Y*) instead of the previous sample {Xi,Yi). 

We suggest to repeat several times (let say A^^ times) this bootstrap procedure, and 
to use the empirical distribution of f*{x) ~ ^{x) fo^' bandwidth selection purpose. 
Precisely, the bootstrapped bandwidth is defined as follows: 

Definition 1 Given Nb replications of the above described bootstrapping scheme, 
and given a fixed set H of bandwidths, the bootstrapped bandwidth h* is defined by: 



Of course, this procedure has still to be validated theoretically (see discussion in 
Sectional), but we will see in the next Sections 15.31 and 15.41 how it behaves both on 
simulated on and real data samples. 

5.3 Some simulations 

The aim of this section is to look at how the automatic bootstrapped bandwidth 
constructed in Definition ^ behaves on simulated samples. We construct random 



h* = h*{x) = argmin 




curves in the following way: 



X{t) = sin{ujt) + (a + 2n)t + b, t e (-1, +1), 



where a and b (respectively uj) are r.r.v. drawn from a uniform distribution on (0, 1) 
(respectively on (0,27r)). Some of these curves are presented in Figure [T] below. 
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Figure 1: A sample of 30 simulated curves 



The real response is simulated according to the following regression relation: 



and where e is drawn from a A/"(0, 2) distribution. 

Our experience is based on the following lines. For each experience, we simulated 
two samples: a sample of size n — 100 on which all the estimates are computed and a 
testing sample of size n = 50 which is used to look at the behaviour of our method. 
Also, for each experience, the number of bootstrap replications was taken to be 
Nb = 100. Other values for J and Nb were also tried without changing the main 
conclusions. To improve the speed of our algorithm, the bandwidth h is assumed to 
belong to some grid in terms of nearest neighbours, that is 



Y = r{X) + e. 



where 




h = h{x) e {hi, /i32} = H, 



(7) 
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where hk is the radius of the ball of center x and containing exactly k among the 
curves data Xi, . . . Xiqq. Concerning the other parameters of our study, the kernel 
function K was chosen to be K{u) = 1 — m^, u G (0, 1) and the norm 1 1. 1 1 was taken 
to be the L2 one between the first order derivatives of the curves. 

We computed, for the 32 different values of h, the average (over the x's belonging 
to the second testing data sample) of the true error (^(x)^^(x))^ a-iid of its bootstrap 
approximation (r*(x) — r(x))^. Finally, this Monte Carlo scheme was repeated 
J = 100 times and the results are reported in Figure |21 (only 25 among the 100 
curves are presented to make the plot clearer). 




5 10 15 20 25 30 5 10 15 20 25 30 



number of neighbours number of neighbours 

Figure 2: Simulations: True (left) and Bootstrap (right) Errors as functions of h 

It appears clearly that both the theoretical quadratic loss and its data-driven 
bootstrapped version have the same convex shape. This convex shape is directly 
linked with the asymptotic expansion given in Theorem ^ before: large values of 
h give high bias, while small values of h lead to high variance. These results are 
quite promising in the sense that the similarity of the shapes of both sets of curves 
presented in Figure |21 let us expect that the bootstrapped bandwidths will be closed 
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from the optimal ones. To check that point , we computed the theroretical minimal 
quadratic loss (that is, the error obtained by using the best bandwidth) and we 
compared it with the error obtained by using the boostrapped bandwidth h* . This 
was done for each among the J = 100 experiences, and the results are reported 
in Figure El which gives mean, variance and density estimates of these two errors. 
Undoubtedly, these results show the good behaviour (at least on this example) of the 
bootstrapping method as an automatic bandwidth selection procedure. Of course, as 
discussed in Section IHl theoretical support for this functional bootstrap bandwidth 
selection rule is still an open question. 




o.o 0.5 1 O 1 .5 2.0 O.O 0.5 1 o r .5 2.0 

"Theoretical Optimal MSE IVISE with data-driven bootstrapped bandwidi 

Figure 3: Simulations: MSE with optimal (left) and bootstrapped (right) bandwidth 



5.4 A real data chemometric application 

Let us now quickly show how our procedure is working on real data. These data 
contain of 215 spectra of light absorbance (Zj, i = 1, . . . 215) as functions of the 
wavelength, and observed on finely chopped pieces of meat. We present in Figure 0] 
the plots of the 215 spectra. 
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wavelenghts 

Figure 4: Spectrometric curves data 

For each spectral cruve corresponds some real response Yi which is the percentage 
of fatness, and our aim is to study the regression relation existing between the real 
variable Y and the functional one Z. These data have been widely studied and, 
inspired by previous studies (see Ferraty & Vieu, 2006) we decide to apply the 
functional kernel methodology to the curves X = Z", and by taking as norm ||.|| 
between curves the usual L2 norm between the second derivatives of the spectra. 
The kernel function K was chosen to be K{u) = 1 — u"^, u G (0, 1). Along our study 
we splitted the data into two subsamples. A first subsample of size n = 165 from 
which our estimates are computed, and a testing sample of size 50 on which they 
are applied. 

In a first attempt, we used the automatic bootstrapping bandwidth selection rule, 
where H was defined as in ((7j). We present in Figure El the shape of the Bootsrapped 
Mean Square Error as a function of the number of neighbours (and thus, as function 
of the bandwidth). 

The same convex form as for the simulated data appears. This form matches the 
theoretical results obtained in Section 13.11 with high bias for large values of h and 
high variance for small bandwidths. These bootstrapped errors are, in this example, 
minimal for the value k = 8. That means that, for each new curve x to be predicted, 
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Figure 5: Spectrometric data: Bootstrapped Errors as function of the bandwidth 

the data-driven bootstrapped bandwidth h*{x) is such that there are exactly 8 
curves-data which are falhng inside of the ball of radius h*{x)- 

These bandwidths lead to completely automatic data-driven fat contents pre- 
diction. For instance, we present in Figure IHl the fat content predictions for the 
50 spectra in our testing sample. In order to highlight the nice behaviour of our 
prediction algorithm, Figure IHl plots the predicted values as functions of the true 
ones. 




10 20 30 40 

True fat content 



Figure 6: Spectrometric data: Predicted values on the testing sample 
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5.5 About implementation of the method 

The implementation of the method was performed by using the Spins routine funo- 
pare. kernel which is included in the package npfda. This package will go with the 
monograph by Ferraty & Vieu (2006). This Spins package (as well as a similar R 
package) and the spectrometric dataset (as well as several other curves datasets) 
will be put in free access on line in the next future. By that time, programs and 
data are available on request. 

6 Conclusions and open problems 

This paper completes the recent advances existing in kernel nonparametric regression 
for functional data, by giving not only the rates of convergence of the estimates but 
also the exact expressions of the constant terms involved in these rates. These theo- 
retical results deal with mean squared errors evaluations and asymptotic normality 
results. As explained in Section El these new results open interesting perspectives 
for applications, such as for instance data-driven automatic bandwidth selection and 
confidence band construction. 

We support the idea that bootstrap methods offers interesting perspectives for 
the functional context. One of them is illustrated by our functional version of the 
wild bootstrap for selecting the smoothing parameter. We have observed nice results 
for bandwidth selection on some simulated and real datasets. It should be pointed 
out that theoretical support for this functional wild bootstrap bandwidth selection 
rule remains an open problem. Our guess is that it should be possible to extend 
to functional variables some results stated in finite dimension (for instance those in 
Hardle & Bowman, 1987), but this has still to be proved. 

Another direct application of our result concerns the construction of confidence 
bands. Once again, the problem of estimating the constants involved in the asymp- 
totic normal distribution can be attacked by the wild bootstrap track described be- 
fore. One possible way for that would be to try to extend standard finite-dimensional 
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knowledge (see for instance Hiirdle & Marron, 1991, or Hardle, Huet & Jolivet, 1995) 
to infinite-dimensional variables. 



7 Appendix: Proofs 

In the following, given some M-valued random variable U, will stand for the 
probability measure induced by U. To make its treatment easier, the kernel estimate 
f will be decomposed as follows: 



dix) 
fix) 



where 



and 



dix) 



1 " /I 



- x\ 
h 



1 " /I 

/M-^E^-(J 



Xk-X\ 
h 



7.1 Proof of Theorem 1 

The proof is split into two parts: computations of the bias and of the variance of 
the estimate. Each part is decomposed in technical lemmas that will be proved in 
Section 17.31 

- Bias term: proof of Let us write the following decomposition. 



A, 



Efix) (Efix)Y (Ehx)'^ 



with 



and 



Ai = E 



A,, = E 



g{x)(f{x)-Ef ix) 



f{x)-Ef{x)] ?{x) 



(9) 



(10) 
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The first step of tlie proof consists in rewritting the first term in right hand 
side of the decomposition (jH)) in the following way: 



Lemma 1 We have: 

Egjx) 
Efix) 

where 



-r{x) = hif' (0)/ + o(/i), 



^ _ /o tK (t) dPll'^-^ll/'^ (t) 
~ K {t) dP\\^-x\\/h (t) ■ 

In a second attempt the next lemma will provide the constant term involved 
in this bias expression and its limit. 



Lemma 2 We have : 

^_ Kil)-j;{sK{s)yT,is)ds _^ Mo 
K{1)- K'{s)This)ds ^ Ml 

Finally, to finish this proof it suffices to prove that both last terms at right 
hand side of (jH)) are neglectible. This is done in next lemma. 



Lemma 3 We have: 



Ai = 0{{nF{h)y^) andA2 = O {{nF (h))-^) . 



So the proof of © is complete. 

- Variance term: proof of . The starting point of the proof is the follow- 
ing decomposition. This decomposition has been obtained in earlier work by 
Collomb (1976) (see also Sarda & Vieu (2000)) in the finite-dimensional case, 
but since the proof is only using analytic arguments about Taylor expansion 
of the function 1/z around 0, it extends obviously to our functional setting: 
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Vara (y) ^9 ix) Cov [g{x)J ix) , 

W) = , / ,2 - 4 r^—rs (11) 

[Ef ix)) [Ef (x)^ 



+ Warf (x) — — -J + o 



Finally, the result (jH)) will follow directly from this decomposition together 
with both following lemmas. 

Lemma 4 We have successively: 

Ef{x)-^K{l)-f K'{s)ro{s)ds = M,, 
Jo 

Eg{x) ^ r (x) {k (1) - j^K' {s) To (s) ds^ = r (x) M,. 
Lemma 5 We have successively: 

{Varg{x)) = {^l + (x)) (1 + o (1)) , 

C^ot; (? (x) , / (X)) = r (x) (1 + o (1)) . 



7.2 Proof of Theorem 2 



The following lemma states a preliminary pointwise limiting distribution result. This 
lemma stems from the bias and variance expressions obtained along Theorem 1; it 
will be proved in the next subsection. 

Lemma 6 We have: 

(r (x) - r (x) - 5n) AT (0, 1) . (12) 
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Because of standard Glivenko-Cantelli type results, we have 

F{h) p 
F{h) ^ ' 

and this is enough, combined with the resuh of Lemma to get the conchision of 
Theorem 2. 

7.3 Proofs of technical lemmas 



Proof of Lemma [2 To calculate 



-r{x) = — 



Ef ix) 



EK 



\x-x\ 



(13) 



note first that 



E 



{Y-r{x))K 



\x-x\ 

h 



E 
E 



V{\\X-X\\)K 



\x-x\ 

h 

x-x\ 
h 



Moreover, it comes: 



E 



V{\\X-X\\)K 



\x-x\ 



V(0) tK{t)dP\\^'''\\^^{t) + o{h) 



the last line coming from the first order Taylor's expansion for around 0. 
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For the denumerator in ()13p we have 



E[K{\\X- x\\ /h)] = j K{t) dpII^-xllM (t) . 



Finally, it appears clearly that the first order bias term is hip' (0) /. 



Proof of Lemma [2J We note that 



tK{t)=K{l)- {sK{s))'ds. 



Applying Fubini's Theorem we get 



tK (t) rfpl'^-^^ll/'^ (t) = K{1)F (h) - (^j^^ {sK {s))' ds^ rfP"^-^"/'^ (t) 

= K{l)F{h)- f {sK{s))' F{hs)ds. 
Jo 

Similarly, we have 

K {t) cipll^->^ll/'^ (t) = K{1)F (h) - f K' (s) F {hs) ds. (14) 
Jo Jo 

So the proof of this lemma is finished by applying the Lebesgue's dominated 
convergence theorem since the denumerator may be easily bounded above by 
K (1) > {K being decreasing). 

Proof of Lemma |1] The first assertion follows directly from ()14|) . while the 
second one can be proved similarly according to the following lines: 



EYK 



\x-x\ 



E 



E{Y\X)K 



\x-x\ 



h 



(r {x) + o{l))E 



K 



\x-x\ 



Proof of Lemma HJ We write the variance of / (x) as: 
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Varf (x) ' 



h 



nF^ (h) 

and note that, as for getting (HD), it holds: 

EK' (^^^) = F (h) (^K' {!)- j'^ {K^y (s) T, (s) ds^ . (15) 
The first assertion of lemma |^ gives 



At last 



Varfix)) -{nF {h))-' {k' (1) - {K')' {s) r, (s) ds^ , (16) 



which finishes the proof of the first assertion of our lemma. 

The same steps can be followed to prove the second assertion. We write 



{Varg{x)) 



nF^ (h) 



The second term at right hand side of this expression is treated directly by 
using the second assertion of Lemma IH while the first one is treated as follows 
by conditioning on X: 



h J ' ' \ h J ' \ h 

The continuity of and of (T^(.) insure that 

y^rgix) = {-I + (x)) EK^ (^^) (1 + o (1)) 
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Combining this result with ()15|) allows to finish the proof of the second asser- 
tion of our lemma. Let us deal now with the covariance term : 



Cov [g{x), f{x) 
1 



nF2 (h) 



EYK^ I /vn \ _;ex { "'" . 1 EYK 



X~x\\\ ^,^f\\X-x\\\ ^,,,,f\\X-x 



h \ h \ h 



The last two terms were computed before, while the first one is treated by 
conditioning on X and using continuity of r: 

EYK' (^^X^) = {r ix) + 0(1)) EK' (^^^) • 

The proof of this lemma is now finished. 

- Proof of Lemma |^ Both assertions of this lemma are direct consequences 
of Lemmas |3] and El 

- Proof of Lemma ini On one hand, (jH)) and Lemma El allows us to get 
r[x) - Er[x) = t^— - + o 



fix) Efix) \V^^, 
On the other hand, the following decomposition holds: 

g{x) Egix) _ ^^ix) - Eg{x)) Efix) + {Efix) - fix)) Eg{x) 
Tix) ~ Efix) ~ fix) Efix) ' 

Using Slutsky's theorem and Theorem (ll we get 



'nFih)Ml 



2 



2M {rix)-Erix))-^MiQA) 



£ 



noting that r (%) — E'rix) can be expressed as an array of independent centered 
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random variables (and the Central Limit Theorem applies). Let us remark that 



r{x)-Er{x) = r{x) - r{x) - B„ 
which achieves the proof of this lemma. 

7.4 Proof of Propositions [T] and [2] 

- Proof of Proposition and ii: Obvious. 

- Proof of Proposition l^-^*^: We have: 



r,(s)=exp(-^fl-l) ]sj. 



If s = 1 we have Th{l) = l,V/i, while if s = we have t/^ (1) = tq (1) = 
F(0) = 0,Wh. To complete this proof it suffices to note that, for s G (0, 1), we 
have: — 1 > 0, and so we have (s) — as /i — ^ 0. 

Proof of Propostion : For any s > and any /i < 1, we have: 

\lnh\ Inh 1 

'^^ ~ \lnh + \ns\ ~ Inh + lns ~ TTW 

I ' Inh 

and so we have to(s) = 1,Vs > 0. To complete this proof it suffices to note 
that Th (1) = To (1) = F(0) = 0, V/i. 

Proof of Propostion [21-i: By simple integration by parts we arrive at 

Mo= [ sK{s)To{s)ds, 
Jo 

where = tq. Because of H2 and because tq is non decreasing, there exists 
some nonempty interval [a, b] C (0, 1) such that both K and Tq do not vanish 
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on [a,b], and therefore we arrive at: 

Mo> [ sK{s)To{s)ds > 0. 



- Proof of Propostion ^i: Obvious. 
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